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^ ■ Abstract 



We calculate the 0{as) corrections to the Z'^-penguin and electroweak box diagrams rel- 
' evant for non-leptonic AF = 1 decays with F = S,B. This calculation provides the 

complete 0{awCis) and O^awCtsSiia^ 9wmf) corrections {a^ = a/sin^Ow) to the Wil- 
son coefficients of the electroweak penguin four quark operators relevant for non-leptonic 
K- and B-decays. We argue that this is the dominant part of the next-next-to-leading 
(NNLO) contributions to these coefficients. Our results allow to reduce considerably the 
uncertainty due to the definition of the top quark mass present in the existing NLO cal- 
culations of non-leptonic decays. The NNLO corrections to the coefficient of the color 
singlet (y — A) ® (y — A) electroweak penguin operator Qg relevant for S-decays are 
generally moderate, amount to a few percent for the choice mt{fit = ^t) and depend only 
weakly on the renormalization scheme. Larger NNLO corrections with substantial scheme 
dependence are found for the coefficients of the remaining electroweak penguin operators 
Qt, Qs and Qiq. In particular, the strong scheme dependence of the NNLO corrections to 
Cg allows to reduce considerably the scheme dependence of Cs{Qs)2 relevant for the ratio 
e'/e. 
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1 Introduction 



Electroweak penguin operators govern rare semi-leptonic decays such as — > vrz/i^, Kl — > 
7r^e~^e~, B -nvv and B^^ /i/i and contribute sometimes in an important manner to 
non-leptonic K- and B- decays. Among the latter one should mention the CP-violating 
ratio e' je in Kj^ —>■ vrvr decays and electroweak contributions to two-body decays like 
B nK, Bs TT°(f) etc. §. 

The effective weak Hamiltonian for these decays has the following generic structure 

^e// = ^EK5KMa(/x)Q.. (1.1) 

Here G^? is the Fermi constant and Qi are the relevant local operators which, in addition to 
electroweak penguin operators, include current- current operators, QCD-penguin operators 
and magnetic penguin operators. The Cabibbo-Kobayashi-Maskawa factors V^^^a^ 
the Wilson coefficients Ci describe the strength with which a given operator enters the 
Hamiltonian. The decay amplitude for a decay of a meson M to a final state F is simply 
given by (F|K//|M). 

The renormalization scale (/i) dependence of C"^(yu) = (Ci(/i),...) is governed by the 
renormalization group equation (RGE) whose solution is given by 



,7^(^; 



Tg exp / dg^ 



asitiw) ^ Pig's, 



Cifiw) , (1.2) 



where Tg is the (^-ordering operator, fiw = C){Myy) and fi is (9(mb) and 0{1 GeV) for 
B-decays and K-decays respectively. P{gs) governs the evolution of the QCD coupling 
constant Qs and 7 is the anomalous dimension matrix which depends on the QED coupling 
constant a in addition to = qI/^t^- In what follows we will work to first order in a. 

Now, the initial conditions C{^y/) are linear combinations of the so called Inami-Lim 
functions [0| such as BqIxi) with Xt = mf/M^ resulting from box diagrams, Colxt) from 
Z°-penguin diagrams, Do{xt) from the photon penguin diagrams, Eo{xt) from the gluon 
penguin diagrams etc. The full list of functions including also those relevant for AS = 2, 
AB = 2 transitions and for radiative B- decays can be found in 0. We will give explicit 
expressions for some of these functions below. 

As shown in |^ any decay amplitude can then be written as a linear combination of 
the Inami-Lim functions to be denoted by Fj°^(xt) 

A(decay) = Po(decay) + P,(decay) ) (xt) , (1.3) 

r 

where the sum runs over all functions contributing to a given decay and Pq summarizes 
contributions from internal up and charm quarks. The process dependent coefficients Pr 
include the effect of the renormalization group evolution from down to fi given in 
( p..2|) as well the matrix elements of the operators Qi. On the other hand the Inami-Lim 
functions Fj^^^Xt) are process independent. That is for instance the functions Bq and Cq 
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enter both the semi-leptonic rare decays K — > -kvv and e' / e which is related to non-leptonic 
decays ixix . 

This Penguin-Box Expansion is very well suited for the study of the extensions of the 
Standard Model (SM) in which new particles are exchanged in the loops. We know already 
that these particles are relatively heavy and consequently they can be integrated out 
together with the weak bosons and the top quark. If there are no new local operators the 
mere change is to modify the functions F^^\xt) which now acquire the dependence on the 
masses of new particles such as charged Higgs bosons and supersymmetric partners. The 
process dependent coefficients Pq and remain unchanged unless new effective operators 
with different Dirac and color structures have to be introduced. 

Now, the universal Inami-Lim functions F^^\xt) result from one-loop box and penguin 
diagrams without QCD corrections and it is of interest to ask how these functions are 
modified when 0{as) corrections to box and penguin diagrams are included. 

The interest in answering this question is as follows: 

• The estimate of the size of QCD corrections to the relevant decay branching ratios. 

• The reduction of various unphysical scale dependences. In the case studied in this 
paper, this is in particular the dependence of QCD corrections on the scale ^it at 
which the running top quark mass is defined. As some of the functions F^^\xt) 
depend strongly on m^, their /x^ dependence may result in the uncertainties as large 
as ±15% in the corresponding branching ratios. Only by calculating QCD corrections 
to box and penguin diagrams can this dependence be reduced. 

• The universality of the top dependent functions can be violated by 0{as) corrections. 
For instance in the case of semi-leptonic FCNC transitions there is no gluon exchange 
in a Z°-penguin diagram parallel to the Z'^-propagator but such an exchange takes 
place in non-leptonic decays in which all external particles are quarks. The same 
applies to box diagrams contributing to i^' ^ vrz/z/ and e' /e respectively. It is of 
interest then to find out whether the breakdown of the universality is substantial. 

• Most importantly, however, the inclusion of 0{as) corrections to penguin and box 
diagrams relevant for non-leptonic decays justifies the simultaneous inclusion of par- 
ticular next-next-to-leading (NNLO) QCD corrections to the renormalization group 
transformation in ( [1. 21 ). In the case of the electroweak penguin operator Qg (see 
( p.5| )), relevant for e' /e, this results in a welcome renormalization scheme depen- 
dence of CsifJ') which in turn allows to reduce considerably the renormalization 
scheme dependence of e'/e present at NLO. 

So far the following 0{as) corrections to box and penguin diagrams have been calcu- 
lated: 

• 0{as) corrections to Z°-penguin function Cq and to the box-diagram function Bq in 
the case of rare semi-leptonic decays like K — > nuu, B ^ jiji etc. |^-^ 
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• 0{as) corrections to chromomagnetic and magnetic penguins relevant for radiative 
decays B X,7 and B X,l+l- 

• 0{as) corrections to the photon penguin function Dq relevant for B — > Xgl^l^ and 

• 0{a1) corrections to the matching conditions of all the operators Qi-e, including 
the 0{as) contributions to the gluon penguin function Eq 

The purpose of the present paper is the calculation of 0{as) corrections to Z°-penguin 
and box-diagrams relevant for non-leptonic decays, such as two-body B-meson decays and 
Kl — s> vrvr {s' /e). This will allow to reduce the /xr dependence in the NLO expressions 
present in the literature, to investigate the breakdown of the universality of the relevant 
Inami-Lim functions and to study the renormalization scheme dependence at the NNLO 
level. 

As we will discuss explicitly in Section 3 the QCD corrections calculated here are a 
part of the complete next-next-to-leading (NNLO) corrections to non-leptonic decays in 
a renormalization group improved perturbation theory. In order to complete the NNLO 
calculations of the relevant Wilson coefficients one would have to calculate 0{a1) correc- 
tions to QCD penguin diagrams and in particular perform three-loop calculations 
0{aa1) of the anomalous dimensions of the full set of operators, which is a formidable task 
and clearly beyond the scope of our paper. However our calculation is sufficient to obtain 
the complete O^an^as) and the 0{aw(y.s sin^ O^iTif) corrections to the Wilson coefficients 
Cj-ioifJ') of the electroweak penguin operators where ayy = a/ sirL^O^. It is also suffi- 
cient to investigate the issue of the /it-dependence of these coefficients and of its reduction 
through 0{as) corrections calculated here. Finally it allows to analyze the breakdown of 
the universality in the Inami-Lim functions related to Z'^-penguin and box diagrams. 

Our paper is organized as follows. In Section 2 we recall the known effective Hamilto- 
nian for A^* = 1 decays at the NLO level. We list the contributing operators and give the 
expressions for the Inami-Lim functions. In Section 3 we discuss our paper in the context 
of a complete NNLO calculation, we motivate our approximations and we outline the strat- 
egy. In Section 4 we elaborate on the renormalization scheme dependence. The calculation 
of the gluonic corrections to the penguin diagrams and to the electroweak box diagrams 
is described in Sections |^ and respectively. In Section 7 we collect the results in terms of 
0{awCXs) contributions to the Wilson coefficients Cj-ioi^M^) and discuss their numerical 
relevance at various scales, as well as the residual scale and scheme dependences. Finally, 
in Section 8 we summarize our paper and we briefly discuss the impact of our flndings on 
the phenomenology of non-leptonic decays. 



2 Notation and Conventions 

In this section we establish our notation and recall some deflnitions that will be useful in 
the rest of the paper. We give explicit formulae for AS = 1 decays. It is straightforward 
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to transform them to the AS = 1 case. The effective Hamiltonian for = 1 transitions 
can be written as 0: 

10 

(2.1) 



G 



HesiAS = 1) = v*v,, Y.c^{^i) 



where we have dropped the terms proportional to V*gVud which are of no concern to us 
here. In Ci{fi) = Vi^jj). The operators Qi are given exphcitly as follows: 
Current— Current : 



Ql = {SaU(^)v~A {Upd^)v-A 

QCD— Penguins : 

Qz = {sd)v-A {^l)v-A Qi 

q=u,d,s 



Q2 = {su)v-A {ud)v-A 



[Sadi3)v-A {'J0'Ja)V-A 
q=u,d,. 



Qb = {sd)v-A Y i<l<l)v+A 
q=u,d,s 

Elect roweak— Penguins : 

3 

= 2 isd)v^A Y ^1 iQl)v+A 
q=u,d,s 



Qq — {Sadi3)v-A Y (?/39'ci)v+A 
q=u,d,s 



V+A 



q=u,d,! 



3 3 

Qg = - {sd)v-A Y (^qi^l)v-A QlO = {Sad(j)v-A Y ^9 {Wo)v-A 
q=u,d,s q=u,d,s 



(2.2) 

(2.3) 
(2.4) 

(2.5) 
(2.6) 



Cio(M, 



11 as{M, 
2 



1 



47r ' 
11 asjM^ 
6 47r 



35 a 
18 4^ 



Here, Cg denotes the electrical quark charges reflecting the electroweak origin of Qy, . . . , Qiq. 
The initial conditions for the Wilson coefficients Cj at /i = obtained from the one-loop 
matching of the full to the effective theory are given in the NDR renormalization scheme 
as follows IT^ : 

C2(M, 
C4(M^ 



■^47^^o(x0 + ^ [25o(xO + Co(xO] 
247r OTT 



Sir 



247r 



Eo{xt) , 



67r 



Eoixt) 



— sm 9^ 

OTT 

0, 



4Co(xi) + Do{xt] 



sin' eU'^Coixt) + Doixt)) + lOBoixt) - ACoixt] 



(2.7) 
(2.8) 
(2.9) 
(2.10) 
(2.11) 

(2.12) 

(2.13) 
(2.14) 
(2.15) 
(2.16) 
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where we have introduced a^y 
We recall that 



(7^/47r = a/ sin^O^y. g is the weak coupling of SU{2)l. 



Bo{xt) 
Co{xt) 
Doixt) 
Eoixt) 
Do{xt) 



1 
4 

Xt 



Xt 



I- Xt 

Xt-Q 

Xt-l 

- In + 



+ 



+ 



Xt In Xt 
{xt - If 
3st + 2 

{Xt - 1)2 



In Xt 



4 , -19a;? + 2hx1 



9 



36(xt - 1) 



+ 



2. Xt{lS - llxt - x1) x?(15 - 16xi + 4x? 



xf{5xf — 2xt — 6) 
I8{xt - 1)4 

„2/ 



Inxj 



-- Inxt + 



12(1 -xt)3 



+ 



6(1 -Xt) 



Eo{xt) = Eo{xt) - - . 



(2.17) 
(2.18) 
(2.19) 
(2.20) 
(2.21) 



BQ{xt) results from the evaluation of the box diagrams, Co{xt) from the Z^-penguin, 
Dolxt) from the photon penguin diagrams and EQ{xt) from QCD penguin diagrams. The 
constants —4/9 and —2/3 in ( p.21| ) are characteristic for the NDR scheme. They are absent 
in the HV scheme. For fi ^ non-vanishing and Cio are generated through QCD 
effects. The formulae ( p.7|) -( ^TT6D apply also to the /S.B = 1 case with the appropriate 
change of fields in Qi. 

Let us next recall that in the leading order (LO) of the renormalization group im- 
proved perturbation theory in which (tts^)"" and atlagt)"" terms with t = In {M^/fi"^) are 
summed only C2{Myy) = 1 is different from zero. In particular C-r^io^M^y) = 0. The 
initial conditions given in ( |2.7| )-( prT6| ) are appropriate to next-to-leading order (NLO) in 
which as{ast)" and a{ast)^ terms are summed. In the next-next-to-leading order (NNLO) 
in which a^(ast)" and aas(ast)" are summed, Ola"^) terms in the initial conditions of 
Ci-q{M^) and 0{aas) terms for all Cj(Miv) have to be included. In the present paper 
we will calculate the dominant 0{aas) corrections to the coefficients C7_io(Myi,.) of the 
electroweak penguin operators. As we will discuss below this will be sufficient to sum the 
dominant contribution of the aas{<yst)"' logarithms. 

Finally we should stress that among 0{aas) terms we distinguish between O^awCts) 
and O^awds sin^ 9^) terms for reasons to be explained in detail below. 



3 General Structure at NNLO and Strategy 

Our aim is to compute the 0{as) corrections to the Z°-penguin diagrams and electroweak 
box diagrams relevant for non-leptonic AF = 1 decays. This calculation constitutes 
only a part of the complete computation of the Wilson coefficients Cj(/i) (i=l,...10) at 
NNLO in the renormalization group improved perturbation theory. On the other hand, 
as we will now demonstrate, our results combined with the known 0{as) and 0{a'^) 
anomalous dimensions of Qi provide the complete 0{anras) corrections to the Wilson 
coefficients C7_io(/i) of the electroweak penguin operators not suppressed by sin^ and 
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the 0{a y^ois sin 6yy) corrections quadratic in to these coefficients. These corrections 
turn out to be by far the dominant contributions to Ci^iolfi) at the NNLO level. 

In order to prove these statements it is instructive to describe the computation of the 
Wilson coefficients Cj(/i) including LO, NLO and NNLO corrections. Generalizing the 
standard procedure at NLO [HIH to include NNLO corrections we proceed as follows. 



Step 1: An amplitude for a properly chosen non-leptonic quark decay is calculated 
perturbatively in the full theory including all sorts of diagrams such as QCD penguin 
diagrams, electroweak penguin diagrams, box diagrams, W-boson exchanges and QCD 
corrections to all these diagrams. The result including LO, NLO and NNLO correction is 
given schematically as follows: 



47r I 47r / ^ Ati An An '''' 



= {QiM^)fC{M^), (3.1) 



where (Q^) is a ten dimensional column vector built out of tree level matrix elements 
of the operators Qi. The superscripts (0), (1) and (2) denote LO, NLO and NNLO 
contributions, respectively. The 0{aas) corrections include O^awC^s), O^awC^sSin^ 0^^) 
and O^a^as sin'^ 9^) terms. 



Step 2: In order to extract the coefficients C{M^y) from one has to calculate 
the matrix elements of Qi between the same external quark states as in Step 1. This in- 
volves generally the computation of the operator insertions into the current-current, gluon 
penguin and photon penguin diagrams of the effective theory (W, top and Z° have been 
integrated out) together with QCD and QED corrections to these insertions. Including 
LO, NLO and NNLO corrections one finds 



(g(M,)) 



An ^ I An ^ An^ ^ An An 



(3.2) 



with fj being lOx 10 matrices. As W and Z° have been integrated out only 0{a^as sin^ 6^ 
terms are present in 0{aas) corrections. 

Step 3: From ( |3.1| ) and (|3.2|) we extract 



c{M^) = c^^) + ^AMAci^) + (^AMA\" cfY^^^^^ (3.3) 

^ ' An \ An An An An 



where 



^(0) ^ ^(0) ^ (0,1,0,...,0)'^, (3.4) 
(7« = i«-f«"(7(°), (3.5) 
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C« = A«-r«'c(°), (3.6) 
(7(2) = (3_7) 

r^(2) _ T(2) _ fXl^f^il) _ f.{2)TA{fi) _ ^(l)T^(l) /o ox 



The electroweak penguin components of C'^'^ and the components of which con- 
tribute to the coefficients of the operators Qi-q are given exphcitly in (|2:^)-(CT). Cf^ 



has been calculated in |13], but it contributes only to 0{a1) corrections to the coefficients 

of Qi_6. 

Among the four terms contributing to C'^J only A^^) and r'^p'^ are of interest to us 
as these are the only ones contributing to 0{awOis) corrections and to 0{awOis sin^ 6^-) 
corrections quadratic in irit which we aim to calculate. The third and fourth term in ( p.8|) 
contribute only to O^awO-s sin^ 6^) corrections. The third term is mj-independent and 
unknown. The last term can be extracted from the known one-loop results but it has no 
dependence quadratic in and will not be included here. 

The main purpose of this paper is then the two-loop calculation of Z^-penguin and box 
diagrams giving A^^) . The contribution from the effective theory can be extracted 

from the known one-loop results. In the case of the Wilson coefficients of electroweak 
penguin operators a simplification occurs as only the operator insertions in current- current 
topologies in contribute. Were we interested also in the coefficients of the QCD- 
penguin operators, also the insertions into QCD-penguin and QED-penguin topologies 
would have to be retained. 

Step 4: We use next the renormalization group transformation to find 

= f/(/i, M^, a)C{M^) , (3.9) 

where 



with the pure QCD evolution given by 



C/(/i, M^) = [/(°)(/^, M^) + U^^\ix, M^) + f/(2)(/i, 



(3.10) 



(3.11) 



The matrices f/'-*-' and -R^*) are functions of the anomalous dimension matrices of the op- 
erators in question and of the QCD beta function. Explicit expressions for U^'^\ tj^^\ 

and can be extracted from [l^, ITsf. R^"^^ and f/*^^-' are not known as they re- 
quire the evaluation of the three-loop anomalous dimension matrices (9(aa^) and 0{al), 
respectively. From the point of view of the expansion in in the renormalization group 
improved perturbation theory, f/(°\ f/^^^ and JJ^'^^ are 0{1), 0{as) and (9(a^), respectively. 

^(1) and ^(2) are C(l/a,), C(l) and 0{as), respectively. 

Inserting ( |3.1CI| ) and ( |3.3| ) into ( p.9| ) and expanding in we find 



7 



a. r 



(3.12) 



Cs{\i) results from 0{(is) and 0{p?^ terms in C{M^^ and the QCD evolution 

U{lJ,, M^). Ci{fi) results from terms 0{a) and 0{aas) in C(M^y) and tf{^,M^y). Fi- 
nally Cii{ii) is found by taking the contributions C(l), 0{as) and in C(M,;i,) and 
performing renormalization group transformation using I&\ Explicitly we have: 



47r 



47r 



^(0) _^ Qi^(^w) ^(i) 



+ C/(2)(//,Mv,)(7(°), (3.13) 



c^s(M^) 
47r 



a?) 



+ t/«(/.,M^)(7«, 



(3.14) 



^(0) ^ ^s{M^) ^ii) ^ ( ois{M, 



cm 



(7^ + 



An 



An 

+ ^(^)(/.,M^)(7(°). 



(3.15) 



Let us now identify the 0{awOis) contributions to electroweak penguin coefficients, 
calculated in subsequent sections, in this full NNLO result. They are fully contained in 
the last two terms in Ci{ii). These two terms can be schematically decomposed as follows: 



Fi{xt, as) + F2{xt, as) sin^ 9^ 



(3.16) 



The 0{a was) corrections are represented by the first term and our calculation provides 
the complete result for Fi{xt, ag). On the other hand, our calculation gives only a partial 
result for F2{xt,as) which is 0{awasSm^ 9^)- The contributions from gluon corrections 
to photon penguin diagrams and 0{a) corrections to QCD penguin diagrams which both 
contribute to F2{xt, a^) are still missing in the case of non-leptonic decays. Similarly, 
some 0{awO£s sin^ 9w) corrections contributing to the Wilson coefficient functions through 

— * 

Ci[{n) are not known. Yet, as we will argue below all these contributions to the Wilson 
coefficients of electroweak penguin operators are expected to be much smaller than the 
0{avyas) contributions calculated by us. Needless to say there are no contributions to 
Cj-io contained in Cs{[j). 

In order to understand better the dominance of 0{awas) over 0{awassiv? 9w) cor- 
rections, let us look at the NLO result, where the issue concerns the dominance of 0{aw) 
corrections over 0{ay^^ sin^ 9^) corrections. Prom the expressions given in the previous sec- 
tion we see that at NLO C^{M-^) is much larger than C'j{M-^): in units of a-^/^T^ we have 
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CgiM^) = -4.46 and C^iM^) = +0.52. More precisely, Cg{M^) = 2.24 sin^ -4.97, 
namely at the electroweak scale Cg is dominated by the second term, unsuppressed by 
sin^ 9w, while the first one accounts for 10% of the total. The dominant term includes the 
box diagrams and the SU{2)l component of the penguin diagrams, all contributing 
0{g'^) terms unsuppressed by siv? 6^^: it can be called the purely weak contribution. 

As we discussed above, the complete QCD corrections to the coefficients of the elec- 
troweak penguin operators involve the computation of the gluonic corrections to the one- 
loop Z*^ and photon penguins and to the electroweak boxes. The two-loop Z^ penguins and 
boxes can be calculated at 0{aas) setting all external momenta to zero. They are entirely 
responsible for the purely weak contribution to Cg{Myy)^ which is largely dominant at the 
one-loop level, as we have just seen. On the other hand, the calculation of the two-loop 
photon penguins is more involved, essentially because the corresponding diagrams lack a 
heavy mass scale like Mz for the Z° penguins. Very recently, the color singlet component 
of this class of diagrams has been computed in a different context [|1^]: we have verified 
(see Section|^) that its 0{awOisS\v? 9w) contribution to C^^g is small compared to the one 
of the boxes and of the Z^ penguin. The color octet component has not yet been calcu- 
lated. One possible strategy therefore consists in computing the gluonic corrections to Z^ 
penguin and electroweak box diagrams exactly at C(avv«s) and in neglecting all correc- 
tions vanishing as sin^ 0, in particular corrections to the photon penguin. These 
purely weak contributions form a gauge- independent subset. 

Before embarking in a complex two-loop calculation, it is also interesting to see how a 
Heavy Top Expansion (HTE), i.e. an expansion in inverse powers of the top quark mass, 
could approximate the complete result. We notice that at the one-loop level the only 
contributions which are quadratic in the top quark mass originate from the Z^ penguin 
diagrams, i.e. from Cq. This feature persists at the two-loop level 0{awOts)'- restricting our 
analysis to these potentially enhanced contributions would simplify significantly our task. 
However, a closer look at the one-loop Wilson coefficients shows that, despite the fact that 
the HTE of -Bq, C*o and Dq converges rapidly, the leading order of the HTE approximates 
well Cj{M^) but not Cg{M^) (it gives -1.64^ instead of -AA5^; this is due to the 
large coefficient in front of Bq in (|2.171 ), which has no quadratic term in mt). It is therefore 
unlikely that the leading order of the HTE provides by itself a good approximation at the 
two-loop level. On the other hand, keeping all terms unsuppressed by sin^ 9^ (i.e. the 
purely weak ones) together with the leading HTE of the rest in the one-loop expressions 
gives 5% and 0.5% accuracy for Cr{Mw) and for Cg{M^v), respectively. 

In summary, we will compute the QCD corrections to Z^ penguin diagrams and to 
electroweak boxes and exclude all the terms proportional to sin^ 9^ which are not quadratic 
in mt. Our approximation provides the complete 0{aw(Xs) corrections to C7_io(Mvy) not 
suppressed by siii^ 9w, as well as the full O^awCtsSiia^ 9wmf) correction. At the one-loop 
level, the combination of these two approximations reproduces very closely the full results. 
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4 Renormalization Scheme Dependence 



Next we would like to elaborate on the renormalization scheme dependence of the Wilson 
coefficients and its cancelation in physical amplitudes. For the purpose of our calculation 
we will need only the transformation f/(/i, M^^, a) including LO and NLO corrections and 
two NNLO terms to be specified below. Indeed as seen in (gll only f/(°) and IJ^^^ enter 
Ci{n). At NLO we have 



U{fi, M^, a) = (/X, M^)W'{M^) , (4.1) 

where 

U^^\fi, Myy) is the LO evolution matrix for which the explicit expression can be found 
in [0,|14|,|15|. Also expressions for Jg, and Jge can be found there. They are functions of 
one-loop and two-loop anomalous dimensions of the operators in question. 

From (|?Tp-(|]Bp we extract 

U^'\^^,M^) = ^J,?7W(/i,M^)-^i^^?7W(/i,M^)J., (4.4) 
R^'^I^^M^) = -^im^){ij,M^) u(o){^^M^.)l (4.5) 

and suppressing the arguments of ?7*^°^ 

- U^'Use + m'UJg-^^^^lU^'Ug. (4.6) 
Now whereas r^^\ f^^\ Jg and J^e depend on the renormalization scheme of operators it 



can be shown that 

r«" + Js, f«" + J,, (4.7) 



are renormalization scheme independent. At NLO it follows that the scheme depen 



dence of tj^^^C^^ in ( 3.13|) is canceled by the second term in IJ^^^ multiplied by C'-"-*. 



Similarly, the scheme dependences of tj'^^'^C'^'^ and R^'^^C'^^ in (7/ and Cu respectively 
are canceled by /\FS^'> C^'^'> in ( 3.15 ), where Aw^^ represents the three last terms in (4.6). 



The remaining scheme dependences reside in At/'^^^C*^^-' in ( p.l3| ) and A/?*^^^^^^-' in ( |3.15| ), 



where this time AR^^^ represents the first three terms ( [4.6| ) and AU^^^ the ffist term in 
([4.4|). One can verify that the scheme dependence of these terms is canceled by the one 
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of the matrix elements To this end the matrix elements in ( |3.2| ) with 0{as) and 

0{a) terms retained and M^^ — ^ should be used. 

Turning to NNLO contributions in ( p.l4| ) let us concentrate on C^"^^ and in particular on 
the scheme dependent term in ( p.8|) f^^^ C^^-* which is taken into account in our calculation. 
Keeping only this term in (7^^-* and adding the last term in ( p.l4|) we obtain 



- f/(°) (f«" + I) Ci'^ + ^Js f/(°) Ci'^ . (4.8) 

An ^ ^ 47r 



That is the scheme dependence of r^^^^ in Cj^l^ has been canceled by Jg in the second term 
in ( [4.4|) . However as C^^^ is renormalization scheme dependent through f^^^ and Jg in the 



second term in ( [4.8|) is scheme dependent, both terms ( [4. 81) remain scheme dependent. In 
order to cancel the scheme dependence of the first term in (|4.8| ) we would have to know 
other terms in ( p.8| ) which as discussed above we do not know. Fortunately at NNLO 
the remaining scheme dependence in ([4.8| ) does not bother us as it is mrindependent 
and the term (r^^''^+ Jg) r^^^^ does not contribute to 0{ayyas) and to the contributions 
v^a^ sin 6^r) quadratic in rrit considered by us. Thus in evaluating the last term in 
( p.l4| ) we will consistently drop the terms originated in photon penguin diagrams, which 
are scheme dependent. Applying this procedure also to the second term of (|4.8|) one can 



easily verify that the remaining scheme dependence of this term residing in Jg is canceled 
by the scheme dependence of {Q{fi)). This procedure has to be properly implemented 
in our calculation of NNLO matching conditions. C^^-* in ( |4.8| ) is then simply given by 
C7-io{Mw) with Do(xt) removed. Needless to say at NLO the full Cj,^^ should be included. 
Now in Section |^ we will perform the renormalization group evolution in order to calculate 
Cj-ioifJ') A* ^w- From the preceding discussion and ( p.l4|) it should be clear that 



this evolution should include the full NLO evolution modified by the following NNLO 
terms: 

i) 0{ayyag) and O {a^dg siii^ 9wTnj) contributions to Ci_io{M^y). 

ii) The contribution 



A [/<») A« - ^ C/<»' M,,) J, 



(4.9) 



representing the last term in (|3.14|) with C^^-* modified as discussed above. It should 



be emphasized that ([4.9|) is not included in the usual NLO calculations as it is 0{aag 
and belongs to the NNLO contributions. 



5 QCD Corrections to the Z^-Penguin Diagrams 

The first part of our analysis is devoted to the QCD corrections to the penguin diagrams 
originating in exchange. It is convenient to separate these corrections according to 
their structure in color space. Let us denote by 1 and T" = A"/2 the N x N matrices in 
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the color space of SU{N). Since the relevant graphs always involve two quark lines, the 
diagrams containing a gluon attached to a single quark line contribute to the color-singlet 
component, characterized by 1 ® 1, while diagrams where the gluon joins two different 
quark lines contribute to the color-octet component, proportional to T" (g) T". 

As far as leading order and color-singlet two-loop diagrams are concerned, the con- 
tributions of the Z°-penguin vertex to the Qi operators can be described in terms of an 
effective sdZ^ vertex 



sdZ 



9" 



with At = V^gVtd and c 



(167r2) c, 

cos6'h,. The coefficient C{x) at 0{as) can be written as 

C(x)=Co(x) + -^Ci(x), 
47r 



(5.1) 



(5.2) 



where Co(a;), introduced in ( p. 18 ), is the relevant Inami-Lim function and Ci(x), which 
was calculated in [0,131, reads 



Ci(x) 



29x + + Ax^ 
3(1 -x)2 

4x -|- 



X 



35x2 _ 



3(1 



3x^ , 20x2 
mx — 



X 



2(1 -x) 



I 4 

- m X 



X 



Here we have used x 



Li2(l - x) 8x 
fit /Ml, and 
Li2(l-x) = 



>9go(x) 
dx 



Inx 



Int 



X > 0. 



(5.3) 



(5.4) 



The scale fit is the renormalization scale of the MS running top quark mass mt{fit)- We 
recall that Co,i depend on the gauge parameter of the ly-field. This dependence is canceled 
at the level of Wilson coefficients by other contributions, to be considered later on. ( p. 18 ) 
and (|5.3| ) actually hold in the 't Hooft-Feynman gauge, = 1. 

In terms of the effective AS = 1 Hamiltonian, the singlet contribution of the Z^- 
penguin can be written as 



n, 



sdZ 
cfT 



Gf 

71 V * 



Co(Xt) + -^Ci(Xt) 

47r 



E [{n 



o 



LL 



eg Olr 



(5.5) 



q=u,d,s,c,b 



where T3 = ±1/2 is the third component of the weak isospin, is the electric charge 
of the quark flavor q and we have introduced the shorthand notation = sin^^^,. The 
four-quark operators are given by 

Oll = S7^(l - 75) d g7''(l - 75) g = {sd)v-A iqq)v-A , (5.6) 
Olr = S7^(l - 75) d qYi^ + I5) q = isd)v^A iqq)v+A ■ (5.7) 

Cg — 1/6, we can rewrite ( |5.5| ) in the basis of the Qi operators 



Using the identity T| 



n 



sdZ 
eff 



Qs + ^slQr -4cl Qg 



{5.i 
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which modifies at 0{as) the Cq contributions to the Wilson coefficient Cj of (|2.7D-( prTBD . 

This completes the analysis of the color-singlet contribution. We now proceed to the 
calculation of the color-octet corrections to the Z°-penguin diagram, which is absent in 
the literature. We first calculate the two-loop diagrams in the full theory. We will then 
compute the renormalization contributions and finally match the renormalized amplitude 
of the full theory with the result of the calculation of the 0{as) corrections to the effective 
theory as explained in Section 3. The relevant two-loop SM diagrams are displayed in 
Fig.|I|. It is important to realize that this is just a small subset of the 0{g'^gg) diagrams, 
which also include, for instance, all the electroweak corrections to the gluon penguin 
diagrams. Fortunately, because of their ffavor structure, most of them do not project on 
Qj-io and are not interesting for our purposes. Only the diagrams involving a or photon 
exchange across the two quark line^, as in Fig. ^ will contribute to the electroweak penguin 
operators, even if they are originated by a gluon penguin vertex (Cf. Fig.|l| (g)). According 
to the strategy elaborated in Section^ we will compute only the Z° exchange diagrams. 

^We recall that the electroweak corrections to the flavor conserving vertex of a gluon penguin diagram 
vanish as a result of Ward identities . 
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Figure 2: Counterterm diagrams for the color-octet component of the 0{as) corrections to 
the penguin vertex. 

All diagrams in Fig.|l] except (c), (d) and (k) are infrared (IR) divergent. We regulate 
these divergences by the use of a common mass m for the internal light quarks and set 
all external momenta to zero (Cf. H^). It is noticeable that our results for the Wilson 
coefficients are unchanged if we set all light quark masses to zero and regulate the IR- 
divergent Feynman integrals keeping a mass parameter only in the denominators. The IR 
divergences are canceled in the matching procedure by the contributions of the effective 
theory. 

The IR divergent graphs have also ultraviolet (UV) divergences which we regulate in 
n = 4 — 2e dimensions using an anticommuting 75. Some of the UV divergences - the 
ones related to the exchange of a pseudo-Goldstone boson - persist after implementation 
of the GIM mechanism. In a calculation of on-shell amplitudes, they would be canceled 
by external leg corrections. However, the IR regularization we have adopted prevents 
the cancelations among the off-diagonal wave function renormalization of the internal 
quarks which are a prerequisite for this procedure. We are then forced to renormalize 
the amplitude at the diagrammatic level, and we do that by zero momentum subtraction 
of the one- loop sub-divergences (the relevant counterterm diagrams are shown in Fig. |]). 
Specifically, writing the quark two-point function for the j ^ i transition as 

S,,(p) = ^P, + mp') ^P, + mp') KP, + m,P«) , (5.9) 

where P^^h = |(1 =F 75) are the left and right-handed projectors and m^j = m, the 
subtraction involves Sfj(O) and Sfj(0). ^fj{0) is 0{rn?) and can be neglected. This 
subtraction procedure removes the spurious IR sensitivity of the diagrams in Fig. ^ (a)-(b) 
and, in the limit m ^ we are considering, implements the correct LSZ conditions on the 
external legs jl^ . The case of the vertex-subdivergences is easier because the diagrams in 
Fig.|l| (e)-(j) are less IR-sensitive. One can therefore neglect all terms proportional to m 
in ( [5. 91 ) and the subtraction involves only ^P^Sfj(O). For a more detailed discussion of 



the renormalization of off-diagonal quark amplitudes, see |l8|-p0 . 
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We have performed two independent calculations, employing a combination of Math- 
EMATICA ETI routines for the various stages of the computation, from the generation of 



the Feynman diagrams |22|, to the Dirac structure simplification and the two-loop 
integration [p|,|25[|. 

After using the unitarity of the CKM matrix, the renormalized two-loop amplitude in 
the full theory, suppressing the external quark fields, can be written as 



(5.10) 



where the sum runs over k = LL, LR, 1, 2 and the spinor structures are given by 



Tll 
Tlr 
Ti 

T2 

Ts 



L(g)L + R(g)L + L(g)R + R(g)R. 



(5.11) 



cr 



j^R ® YL + ® l^R 



with i?, L = 1±75. In order to project the renormalized amplitudes on the different spinor 
structures, we use the method adopted for example in and reduce the problem to the 
calculation of traces of strings of Dirac matrices. As the amplitudes are now finite, this can 
be done in four dimensions. The coefficients "W^^^ ( furthermore be decomposed 

according to 



>viS(x*,x.) 
>Vf)(xi,x,) 
>Vi')(xi,x,) 



Z{xt,x^) + 6Co{xt) Inxq 



Z(xt,x^ + 6Co(xt) Inx, 



(5.12) 



-(3 + (r|-2e,4)Co(xi) 
(3 + OT|Co(x,), 



where q is the flavor of the lower quark line (Cf. Fig.|I|), Xz = M'^/M^ and the logs of 
Xq = m^/M^ indicate the IR divergences. In analogy to the case described in |T^, the 
structures Ti, T2 and T3 in ( |5.1CI|) are artefacts of the IR regularization procedure and we 
will verify in a moment that they drop out in the matching with the effective theory. For 
instance, if we consistently set the common quark mass to zero in the numerator of the 
quark propagator, W/2'' vanish. In ( ^.12| ) we have left the gluon gauge C, arbitrary and 
set = ^- We have also checked that the Z°-field gauge dependence of the individual 
diagrams cancels in their sum. 

For what concerns the effective theory side, we need the octet-part of the one-loop 
matrix elements of the renormalized operators Oll and Olr in QCD. The calculation 
is performed following the same regularization procedure used for the two-loop diagrams 
and it involves the one-loop diagrams depicted in Fig.^. In principle also insertions in 
the penguin diagrams should be considered. However at the level of the approximations 
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Figure 3: Diagrams contributing to the matrix element of current-current operators at 0{as). 



outlined in Section 3 they do not contribute to the Wilson coefficients of Qr-io- After 
renormalization in the NDR scheme, we obtain 



Hoop 
Hoop 



{OLL)tree + ^T.\CFl(E)iA^^\fi)+T^^T'^Af\fi)] T,, (5.13) 
(OLH)tree + ^E[^^i®l4'^(/^)+^"®™f (5.14) 



with k = LL, LR, 1, 2, 3. The coefficients Ak and Bk which do not vanish are given by 

= -3-2einx, + 2einx^, 

AilifJ') = — 5 + 6 InXg — 6 Inx^ , 

4'^(/i) = -Af\f,) = -2A^\fi) = -(3 + 0, 

^?^(/^) = -3-2einx, + 2einx^, 
^LFiif^) — — 6 In + 6 In , 

Sf^(/i) = Bf\^^) = -2S«(/i) = -(3 + 0- 



(5.15) 



(5.16) 



m. 



The results for Ak can be also obtained from [0 , after taking the limit = rrid 

Unlike the full theory results, ( |5.15|) -( ^.16|) are scheme dependent. For instance the 
constant terms depend on the way 75 is defined in n dimensions — in our case they are 
specific to the NDR scheme. The scheme dependence is generated in the calculation of 
the matrix elements in the effective theory: for example, in the Dimensional Reduction 
(DRED) scheme ||2^,^ there is no constant part in A^^ and S^.*^ but only logarithms. In 
the 't Hooft-Veltman (HV) scheme p9[, the constants in A^lf\jj) are (1, —1) instead of 
(—3, —5), respectively, and in (/i) they are (1, 5) instead of (—3, —7) (see also (3.9)- 
(3.10) of ^B|). ( p.l5| )-( pTT^ ) also depend on the definition of the evanescent operators 
PPj , . It is crucial that this definition follows the one adopted in the calculation of 
the two-loop 0{af) anomalous dimension matrix |15|,|2^. In practice, in our case this 
means that we have to perform the projection on by taking traces in n dimensions with 
anticommuting 75. As a consequence of our choice of IR regularization, and in contrast 
to [|, 32|, this is the only occurrence of the evanescent operators in our calculation. 

We have now all ingredients needed to match full and effective theories. Taking into ac- 
count ( p. 101 ), ( p. 131 ) and ( p.l4| ), we see that the color-octet part of the effective Hamiltonian 
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can be written as 

V q=u,d,s,c,b 

where the coefficients G and H are given by 

G = Z{xt,x^,Sw) + 5Coixt) + 6Coixt)\nx^^ , (5.18) 
H = Z{xt,Xz,Sw) - 7Go{xt) + 6Go{xt)\nx^^ (5.19) 

and we have introduced x^^ = fi'^/M^. The unphysical and ,^-dependent terms obtained 
from the two-loop calculation of (|5.10|) have been canceled by analogous terms from the 
effective theory. Note that the scale fi^ in (|5.18| )- (|5.19| ) is the scale at which the matching 



is performed. The scale /i^^ is not related to the top quark mass renormalization scale /ij 
appearing in ( |5.3| ) although they can be set equal. We will, however, keep them distinct 
in the following. Taking advantage of the identity T| = Cg — 1/6 and using 

T"(g)T" = ^ (5.20) 



2 V N' 



where 1 indicates twisted color indices - like in of (|2.3| ) - we can rewrite (|5.17|) in terms 
of the Qi operators: 

Hf/ = ^^^^t[G{Q3-3Q,)-4 slH {Q, - 3 Qg) - 4 cl G{Q,-3 Q,,] 

(5.21) 

From here one can read the contributions of this class of diagrams to the various Wilson 
coefficients at the matching scale calculated for = 1 in the NDR scheme. 

As we have seen above, at the one-loop level and in the case of the color-singlet 0{as) 
corrections the dependence on ?av? 6w drops out in the functions Co,i(xi). This is a 
consequence of the Ward identity which ensures that the photon exchange diagram has 
no pole, and it is guaranteed because the momentum carried by the Z*^ boson is 
vanishingly small. In the case of the octet contributions the Ward identity does not hold, 
because the momentum carried by the is not small. Indeed, we verify that the sin^ Oy/- 
dependence is not removed from the Wilson coefficient and that the function Z can be 
decomposed into 

Z{xt^ Xz, Svv-) = Zq^Xi, Xz) + Zi(xt, Xz)- (5.22) 



The coefficients Z^ i are complicated functions of Xt and Xz and are given in (|A1|) and ([A2|) 
of the Appendix. As rrit and are now accurately determined, Zq^i can be linearized in 
the vicinity of their central values. Using the latest experimental results mt{mt) = 166 ±5 
GeV, = 80.394 ± 0.042 GeV, and = 91.1867 GeV |3[, we find 



Zo = +5.1795 + 0.038 (mt- 166) + 0.015 (M^^- 80.394), 

Zi = -2.1095 + 0.0067 (mt- 166) + 0.026 (M^- 80.394) (5.23) 
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which reproduce the analytic expressions to great accuracy, better than 0.1%, within 2cr 
from the central values. 

It is also interesting to see how the HTE approximates these two functions. In this 
respect we stress that, although G and H are ^vi'-dependent quantities, their leading HTE 
term is gauge-independent, as Z°-penguins are the only source of contributions quadratic 
in rrit. The contributions quadratic in irit are 

Z^^^ = \xt{l-\nx^) , Zr^ = (5.24) 

so that we find Z^^^ = 2.39 and Zf^^ = at leading order in the HTE. At next-to- 
leading order in the HTE the approximation improves substantially, as we get Zq = 4.62 
and Zi = —1.46, relatively close to the central values of ( |5.23| ). Finally, we notice that the 
leading term of the HTE can be obtained considering only the diagrams involving Yukawa 
couplings of the top quark, as we have explicitly verified. 

In summary, in this section we have calculated the gluonic corrections to the Z°-penguin 
diagrams. The main results are (|5.8|) and (^.21|) . 



6 QCD Corrections to the Electroweak Box Diagrams 



The second part of our analysis concerns the electroweak box diagrams. Again, we will 
consider for definiteness the case of AS" = 1 transitions. Although some 0{as) results 
are available in the literature for the case in which all quark involved in the transition 
are down quarks [|17, 34| and for the case of semi-leptonic transitions the elec- 

troweak box diagrams involving both down and up quark lines require a new calculation 
that we describe in this section. Indeed, it is a fortuitous coincidence that at the one-loop 
level quark box diagrams containing either up or down quarks are described by the single 
function Bo{xt) introduced in ( |2.17| ). As a by-product of this computation we will also be 
able to reproduce all the two-loop box results of [p|,p!7|,p2|,^. 

First, we need to recall some one- loop results necessary for the subsequent discussion. 
The one-loop amplitude for s + d q + q with q = u,c can be written as 



Hoop 



9^ 



Y,X^hjS^^\x.„x,)OLL, 



(6.1) 



Vgjp, Xj = V*gVjd and Oll has been defined in (p^). The function S^'^\xi, Xj) 



where A) 

describes a generic AS" = 1 box with external up quarks and arbitrary internal quark 
masses rriij. Expanding it up to 0{e), it reads 



5(")(x„x,) = 5i"^(a;„x,) + e (x„ x„ x^) + 0{e^) 



(6.2) 



with 



Sq ^ (xj, Xj^ 



16 — 7xiXj 



lQ{xi — l){xj — 1) 16(xj — iy{xi — Xj) 



xf{16 + Xj 



Xi 



^)) 



In Xi + ixi 



, (6.3) 
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40 — ISxiXj 



32{xi - l){xj - 1) 



x]{40 + Xj{3xi - 16)) 



321' 



In Xj 



xfilQ + - 8)) , 2 / 
" Wxi + ixi 



32 ( 



+ 5q Xj) InXp 



(6.4) 



and = m^j/M^. iS*^"^ depends on the ly-field gauge and the above expressions hold in 
the 't Hooft-Feynman gauge. Setting all light quark masses to zero and using the unitarity 
of the CKM matrix, the only relevant combination in the limit e — > is 



5(")(xi,0)-5(")(0,0) = -4i?o(a;i) 



(6.5) 



where Bo{xt) is the Inami-Lim function of ( p.l7|) . Taking advantage of the identity 
J2q=u,c ^LL = iQs + fQg, we obtain the effective Hamiltonian induced by the box di- 
agrams with isospin T3 = 1/2 (up) quarks: 



KfT(T3 = 1/2) = -%^A,i?o(xO (Qs + 2Q,) 



V2 3 



7r 



(6.6) 



The case of the down-quark box diagrams is slightly more complicated in that there 
is a mismatch in the CKM factor between the q = d,s and q = b cases. This implies the 
introduction of two additional operators Qu^u 

Qu = {sd)v-A {bb)v-A , Q12 = {sb)v-A {bd)v-A ■ (6.7) 
Calling S^'^^Xi, Xj) the box function for the down-quark box diagrams, we find after GIM 



{d,s) 
Hoop 



9 

167r2 2MI 



{Xt [S^''\xt,0)-S^^\0,0)]}OLL 



(6.8) 



where we have dropped a term suppressed by A^'^''^'*. On the other hand, in the case of b 
quarks A^'''' is not a suppression factor and we obtain in this case 



m: 



Hoop 



9' 



167r2 AMI 



{\t [S'^'\xt,0)-S^'\0,0)] Oll 

+ Xi'ht [S^''\xt,xt) -2S^''\xt,0) + S^''\0,0)] (Qn + Q12)} , (6.9) 



The function S^'^\xi,Xj) undergoes the same decomposition of (|6.2|) . In the 't Hooft- 
Feynman gauge the coefficients take the form 



<Sq (xi, X 



4 '~l Ob 'i Ob j 



«5 -^J, 



lQ{xi - l){xj 



+ 



xf(4 + Xj(x, -8)) 
16(xi - iy{xi - Xj) 



Inxj + ixi 



(6.10) 



4 + 13x iXj 



32{xi - l){xj - 1) 
x^{4 + Xj{x,-8)) 
32{xi - l)2(xj - Xj] 



+ 



X, 



-A + Xj{3xi - 16)) 



32(a;,; - l)'^{xi - Xi 



In Xi + ixi 



Xj) 



In Xi 



+ Sf^{xi,Xj)hix^. (6.11) 
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Figure 4: Feynman graphs contributing to the 0{as) corrections to the electroweak box 
diagrams. Mirror diagrams are not displayed. 



In n = 4 dimensions the combinations present in (|6.8| ) and (|6.9|) reduce to 

S^''\xt,0)-S^''\0,0) = Boixt) 

and 



S^'\xt,xt)-2S^''\xt,0)+S^'\0,0) ^ ^Soixt) = ^ 



Axt — llx^ + x'l Sx^lnXi 



4{l-xt) 



2(l-xt) 



(6.12) 



(6.13) 



where So{xt) is the box function characteristic of AF = 2 transitions. We see from ( |6.5|) 
and ( |6.12| ) that T3 = 1/2 and T3 = —1/2 box diagrams involve the same function Bo{xt). 
This is true only in n = 4 dimensions and for the 't Hooft-Feynman gauge. We will see in 
the following that there is no such relation at 0{as)- 

Using the identity J2q=d,s,b^LL = ^{Qs — Q9), we can write the contribution to the 
effective Hamiltonian as 

= -1/2) = ^^A,|^i?o(x,) (Q3 - Q9) + ^Af) So{x,) {Qu + Q12) }• (6.14) 
The role of the operators Qn.12 in the RGE evolution of the Wilson coefficients between 



and rrif, has been analyzed in [35]. In the case of e'/e, for instance, they can be safely 
neglected. Their 0{as) corrections are likely to be irrelevant and will not be considered 
in the following. 

We are now in the position to present the calculation of the gluonic corrections to 
the one-loop electroweak box diagrams. The relevant diagrams for the case of isospin 
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T3 = 1/2 are shown in Fig.^. Both color-singlet (c), (d), (g) and color-octet (a), (b), 
(e), (f) diagrams are present. The calculation proceeds along the same lines as the one of 
Section^. Diagrams (e), (f) and (g) present IR divergences which are regulated in the way 
described in the previous section. The origin and the treatment of the UV divergences, 
however, is fundamentally different: diagrams (c), (d) have subdivergences related to the 
quark-gluon interactions and are renormalized in the MS scheme. In fact, it is sufficient to 
renormalize the internal quark masses and to implement the wave function renormalization 
of the external legs. Of course, in the counterterm diagrams the 0{e) parts of ( |6.4| ) and 
( |6.11| ) have to be retained. The renormalized amplitude for the process s + d ^ q + q with 



q = u,c can be written as 



= (iJ)2 4M2^ E ^1'^ A. E [C^ i ® i Vi"' n T« ® r« Vi"' T, (6.15) 

with k = LL, 1,2,3. is the number of colors and Cp = {N"^ — 1)/2N. The spinor 
structures have been introduced in (|5.11|) . Setting = 1, keeping the gluon gauge 



parameter ^ arbitrary, and without making assumptions on the masses of the internal 
quarks, we find 

vtL'\x„x,) = C^^''\x,,x,)-2^lnx,Sl>''\x„x,) (6.16) 

-|-2 ^Inx^ iSg \xi, Xj) + 6 Inx^ i^^'dx- ~^ ftr~) -^j) ' 

Vti:'\x,,x,) = C^^^^\x,,x,)+6lnx,St\x,,x,), (6.17) 
V/"'')(x„x,) = -Vt'\x„x,) = -2Vt''\xi,x,) = -{3 + OSt\xi,Xj). (6.18) 



All remaining V]." ' vanish. The two terms in the second line of ( |6.16D describe the 



scale dependence introduced by the MS renormalization of the external fields and of the 
internal masses, respectively. The functions are independent of the gluon gauge. 

The complete expressions are quite long and can be found in |^ . 



Concerning the effective theory, only the insertion of the operator Oll is relevant in 
this case and the results can be found in the previous section. It is easy to verify that 
all the unphysical spinor structures and the gauge- dependent terms of ( |6.16| )-( |6rT^ ) cancel 
in the matching and we are left only with contributions proportional to T^^. Using the 
unitarity of the CKM matrix and ( p. 201 ), the matching of full and effective theory leads to 
the following contribution to the effective Hamiltonian 

A Hesin = 1/2) = ^^^X, iBUxt) (Q3 + 2 Qg) + B^xt) {Qa + 2 Qio)! . (6.19) 
V2 67r 47r L J 

The functions B'^{xt) and B^{xt) are given by 

B^x^) ^ _2..(23 + 9..) _ 16..(1 - 5x.) _ ..(9 + 23..) 

' Sixt-iy 3(xt-l)3 2(xt-l)3 ^ ' 

rLi2(l - xt) - —Bo{xt) + 'iBo{xt)lnXf,^^ - 32xt — ^ lnx^„ 
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-10 5o(xO-12fio(a;01nx^H'- (6-21) 

Again, these results hold in the 't Hooft-Feynman gauge = 1 and are specific to the 
NDR scheme. As we will discuss in more detail later, the scheme dependence resides in 
the coefficients —38/3 and —10 that multiply Bq^Xi). Notice also that -B" depends on 
both and /ii^. 

We now consider the case of T3 = —1/2 box diagrams. The relevant two- loop diagrams 
are the analogue of the ones shown in Fig.^, although in this case one should also consider 
the Fierz rotated diagrams, which just lead to an overall factor of 2. The renormalized 
amplitude can therefore be written in the same way as in ( |6.15| ), but it is characterized by 
new coefficients V^*^'*^. These coefficients agree with the expressions given in the Appendix 
of [0. After the matching with the effective theory and the implementation of the GIM 
mechanism, we can express the contribution to the effective Hamiltonian of the weak 
isospin T3 = —1/2 box diagrams as 

AKff(T3 = -1/2) = \Bi{xt) {Q, - Q,) + 5f(x,) {Q, - Q,,)] , (6.22) 

y2 OTT 477 L J 

where the functions Bf and Bf{xt) are given by 

, 8-183xt + 47x2 8 + 27a;t + 93x2 xt(27 + 71xt - 2x?) ^ ^ 

= 24(x, - ir 24(x,-l)3 ^""^+ 24(0:, -1)3 

2-3xt-9xJ + xf^. , 2 + x. 19 , 

* *.Li2(l-xJ + --^C(2) + — 5o(x,) 



6xt{xt — ly 6xt 6 

-Bo{xt) Inx^^ + ^^i—^^ l^^Mt ' (6-23) 
Bd/ N 8-23xi S-hxt xt{?> + 2xt) 2 2 + a;, 

BAXt) = ; r ; rTrlnX/H ; TTT In Xf C{2) 

' 8(x,-l) 8(x,-l)2 8(x,-l)2 2xt ^ 



+ o ^1 Li2(l - X,) + ^Boixt) + 35o(x,) Inx^^ . (6.24) 

2xt{xt -ly 2 



As before, the previous expressions are specific to the NDR scheme and are valid for 

Summarizing, in this section we have calculated the NNLO contributions originated in 
electroweak box diagrams. The main results are reported in (|6.19|) and ([6.221). 



7 Numerical Results 

In this section we summarize our results in terms of 0{aw(ys) contributions to the Wilson 
coefficients of the electroweak penguin operators Q7-10 and study their numerical rele- 
vance, both at the electroweak scale and at typical hadronic scales in the NDR and HV 
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schemes. We discuss the reduction of the /it-dependence in the Wilson coefficients and 
the issue of the renormahzation scheme dependence. We conclude with a discussion of the 
universality of the functions X and Y of the Penguin-Box Expansion. 



7.1 Results for the Wilson Coefficients 

Let us collect the results of ( ^^211) , ( |6l9D and Using 



Ho{xt, Xz 



2o{xt,x,) + 5Co{xt) + 6Co{xt)lnXf,^ , 
2:o{xt,x,) - 7Co{xt) + 6Co{xt) In 



(7.1) 
(7.2) 



we obtain the following 0{anras) corrections to the Wilson coefficients of the electroweak 
penguin operators 



aw Ois 
Gn An 

aw as 



2 \ 2 

Sw ( 4Ci(xt) + -HQ^Xt, x^) \ + -s^ Zi{xt, x^ 



^"w Ho{Xt, Xz) + Zi(Xt, Xz 



B^ixt) - Bf{xt) + 2Ci{xt) - -Goixt, 



(7.3) 
(7.4) 

(7.5) 



-s^ (2Ci{xt) - ^^0(0;*, Xz) + ^2i(a;j, x^)^ + ^ s%,Zx{xu Xz 
B'^{xt) - Bf{xt) + Go{xt,Xz) - s'l,{GQ{xt,Xz) - Zi{xt,Xz)) 



'^w Zi{Xt, Xz 



(7.6) 



As we have seen above, there are also contributions to C3.4 which can be extracted from 
(p.8|), ( |5.21| ), ( |6.19| ) and ( |6.22| ). However, any electroweak correction to a gluon penguin 
diagram would contribute at the same order. The subset of diagrams we have computed 
is insufficient for these coefficients. We have organized the results in ( |7.3| )- ([T^ ) according 
to powers of Sw^ It should be clear by now that the zeroth order coefficient is complete 
and gauge-invariant. The same applies to the coefficient of s^, as the only missing part of 
our calculation — the QCD corrections to the photon penguin diagrams — is of 0{awSw) 
and cannot contribute to it. On the other hand, only the leading term of the HTE of the 
coefficient is complete and gauge-invariant. 

As a ffist check of our results, we can verify that the dependence of the NLO coefficients 
Ci on the matching scale fiw and on the top mass renormahzation scale fit is removed by 

^Of course, the argument Xz — {I ~ s?i )^^ of the functions Zq and Zi should also be expanded in 
powers of Sw However, in our approximation the whole ©(m^) term of Zq of (5.24) has to be included 



and can therefore be absorbed in the first term of the expansion. Once this is done, expanding Xz 
becomes numerically irrelevant. 
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NLOndr 






NNLO^Sr 


NLOhv 




C7(M,.) 


0.135 


0.115 


0.116 


0.114 


0.158 


0.142 


C8(M,.) 





0.001 


0.002 


0.002 





-0.005 




-1.091 


-1.004 


-1.002 


-1.014 


-1.067 


-0.963 







-0.019 


-0.024 


-0.019 





-0.003 



Table 1: Wilson coefficients of the electroweak penguin operators at the scale in units a 
in the NDR scheme and HV schemes (see text). 



the scale dependence of the calculated NNLO corrections, up to terms O^an^s^) originated 
by the missing photon-penguins. Indeed, it is straightforward to see that 



X 



(CW(M^) + ^CifiM^)'^ = 0{al) (7.7) 



and similarly for the dependence. This follows from 

^ (7f ) (M^) = -7o- ^ X, (M^) (7.^ 

and 



dx^^ " 2 47r 



Here 7^ = 8 and 7^°-' are the LO anomalous dimension of the top mass and the LO 
anomalous dimension matrix of the operators Qi. Additional dependent contributions 
of 0{aw(y.ss1y) come from the QED induced mixing between the gluon and electroweak pen- 
guin operators. 

The numerical values of the Wilson coefficients at the electroweak scale are reported 
in Table 1, where we compare the NLO and NNLO results. In all numerical calculations 



we employ = 80.394 GeV, = 91.1867 GeV, and OsiM^-) = 0.121 ||. In Table 1 
we furthermore fix fi^ = /it = M^, and consequently adopt mt{Mw) = 175.5 GeV, which 
follows from the experimental value of the pole top mass, rrit = 174.3 ±5.1 GeV. For the 
electroweak mixing angle we use = sin^ 6'^(M2) ~ 0.23145 [^,^. We give three 
different values for the NNLO coefficients in the NDR scheme: NNLO*-^^ corresponds to the 
expressions given in (fO| ) - (|7.6|) , which, as mentioned above, contain some gauge- dependent 
terms calculated in the = 1 gauge. In NNLO*-^-*, instead, we expand the coefficients 
of (7J)-(7^) in inverse powers of rrit and retain only the leading HTE component. To this 



end we recall that 

'4 



Cf™(xt) =^t[^- C(2) + Inx^, - Ina^ij . (7.10) 

The formulation NNLO^^-* is strictly gauge-independent. The QCD corrections modify 
Cj^glMyy) by about —15% and +8%, respectively. The difference between NNLO^^-* and 
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NNLO*^^-* is very small, which is consistent with our expectations about the contributions 
of the QCD corrected photon penguin diagrams. The inclusion of a recent result for 
the color-singlet photon penguin contribution |13] would change the results for CrlMyy) 
very little (by about +3%) and marginally (—0.4%) for Cg{My^,). Following Section^ we 
will not include it here. Finally, the fourth column of Table 1 gives C7_io(M,y) in the 
NDR scheme for the case in which all 0{as) corrections are calculated at leading order 
in the HTE, NNLO^™. The agreement with the third column is relatively good also 
in this case. We also observe that in the C,w = ^ gauge and at fit = the dominant 
NNLO contribution to ACy-io is provided by Ci(xt), the color-singlet corrections to the 
Z^-penguin diagrams. 

The QCD corrections of (|7.3|) - (|7.6|) are specific to the NDR scheme. Using the results 
given in Section^, it is not difficult to find the expressions for the Wilson coefficients of 
([7.3|)-([r6D in the HV scheme: the cofactors of Co(xt) in ([7.1[), ( |7.2| ) become (1,5) instead 
of (5,-7) in NDR; the cofactors of Bq in B'l smd B"^ in ( lOOl) , ( |OTD are (6,-2) instead 
of (-38/3, -10); the cofactors of in Bf and Bf in (lOSi) , ([Ol become (-3/2, 1/2) 
instead of (19/6,5/2). The numerical values of Ct-iq{Mw) in the HV scheme at NLO and 
NNLO for fit = fJ'w = are given in the last two columns of Table 1. Also in this scheme 
the QCD corrections to Cj^q^M^) are 0(10%). We note finally that at NNLO C8^io(M^) 
become non-zero, but still are very small. 



7.2 Reduction of the /it-dependence 

It is interesting to compare the Ht dependence of the Wilson coefficients before and after 
the inclusion of the 0{as) corrections. This is done in Figs. |^ for Cj{M^) and C<j{My^,), 
where we have used the leading log expression for the running mass of the top 

12 
23 

(7.11) 

and employed the NNLO*-^^ expressions in the NDR scheme. 

Despite the fact that the NNLO corrections have not been computed completely, the 
reduction of the scale dependence is remarkable, and is again consistent with the idea that 
the contributions we have calculated are the dominant ones. We also observe that the QCD 
corrections to C7(M^-) and Cq^M^) are particularly small for /ij ^ mt. This has also been 
found in the case of rare semi-leptonic decays As we will discuss below this pattern 
does not apply to Cg and Cio- However, in view of the fact that C%{My^.) = Cioi^M^) = 
at NLO, we will study their /ij-dependence for /i -C Myy. 

In practical applications it is often useful to have simple and compact formulas for the 
Wilson coefficients at the weak scale. For fi^ = and fit = iTit, the NDR coefficients in 
the NNLO'-^^ formulation and in units of a can be written as 

CriM^) = 0.02185 a;,^-^^^2 , Cs{M^) = 0.000718 a;"'^*^^ , 

Cg{M^) = -0.438 , Cio(M^-) = -0.004224 Xj^-^^^^ , (7.12) 



vitifit) = mt{mt) 



as{mt) 
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Figure 5: fit dependence of Ci{M„) and C9(Mjy) at NLO and after the inclusion of the 
NNLO corrections for = Myy in the NDR scheme. 





NLOndr 


NNLOndr 


NLOhv 


NNLOhv 




-0.002 


-0.011 


-0.002 


-0.010 




0.055 


0.060 


0.061 


0.050 




-1.336 


-1.218 


-1.336 


-1.243 




0.277 


0.209 


0.280 


0.260 




-0.030 


-0.032 


-0.028 


-0.037 




0.142 


0.160 


0.151 


0.135 




-1.538 


-1.375 


-1.538 


-1.445 




0.582 


0.441 


0.589 


0.553 



Table 2: Wilson coefficients of the electroweak penguin operators at typical hadronic scales 
fib = 4:A GeV and fix = 1-3 GeV in units of a in the NDR scheme and HV schemes for 

l^^w = lJ^t = Myy (see text). 



which have to be compared with the NLO expressions 

C^'^^iM^) = 0.02268 , C^'''^{M^) = -0.434 , Cl^^^iM^) = . (7.13) 

These expressions reproduce the results of the complete formulas with an accuracy of 0.2% 
or better within two sigmas of the present rrit value. 



7.3 RGE Evolution and Scheme Dependence 

Let us now study the evolution of the coefficients down to a typical hadronic scale. The 
inclusion of NNLO contributions proceeds as explained in Section 3. We will consider two 
cases: the one of B meson decays, for which we will use fib = 4:A GeV and the one of K 
meson decay, corresponding to fix = 1-3 GeV, as used in the analysis of e'/e [p3]. The 
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Figure 6: fit dependence of Csifix) and Cio{fiK) at NLO and after the inclusion of the 
NNLO corrections for /i^^ = My^ in the NDR scheme. 



results for NDR and HV schemes are shown in Table 2 for /it = = M^. We recall 
that part (but not all) of the scheme dependence of Cj^iq is canceled in the evolution 
against similar terms in the anomalous dimension matrix (see e.g. [Q]). As demonstrated 
in detail in Section 4 the renormalization scheme dependence of Cj_io{Mn') discussed 
above is canceled by the first term in ( [4.9|) stemming from the renormalization group 
transformation. The coefficients C{n) are however scheme dependent through the scheme 
dependence at the lower end of the RGE evolution represented by the second term in ( [4.91 ). 
The complete cancelation of the scheme dependence of physical amplitudes occurs only 
with the inclusion of the matrix elements of the operators Qi. We employ ^^(M^) = 0.119, 
corresponding approximately to A*^^) = 340 MeV. In Table |] the entries labeled by NLO 
refer to the NLO case described in Section 2. The entries identified by NNLO, instead, 
correspond to our approximation of the full NNLO result in the form NNLO^^-*. At /i = 4.4 
GeV the shifts due to the new contributions in the NDR scheme and for fit = are 
about +9% for Cg, +9% for Cg, —24% for Ciq. Cj remains very small. At fix = 1-3 GeV 
the situation is similar, although the shifts are naturally more pronounced. In the case 
of the HV scheme the NNLO corrections to Cg are somewhat smaller than in the NDR 
scheme. They are comparable for C7{fj,b) and somewhat larger for Ci^hk)- The strongest 
scheme dependence is observed in the case of Cs and Cio, which is not surprising as Qs 
and Qio are color non-singlet operators. Whereas Cg is enhanced in the NDR scheme, 
it is suppressed in the HV scheme. Cio is suppressed in both schemes but the effect is 
substantial in the NDR scheme and rather small in the HV scheme. 

As explained in Section |^, there are other NNLO contributions that we have ne- 
glected. Some of them are not known, but we can check the magnitude of the neglected 
0{an^assl,) effects from the term as{Mw) / {4:'K)l6^\fi, M^)&^^ in (|3.15| ). It turns out 



that these effects are much smaller than the NNLO contributions we have considered and 
are completely negligible. We also notice that, among the NNLO contributions in ( |3.14| ), 
the one proportional to C'^^ is by far the dominant in the calculation of Cj and Cg. 

Fig. I shows the /xj dependence of C^{iik) and Cio(/ix) at NLO and NNLO order in 
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NLOndr 


NNLOndr 


NLOhv 


NNLOhv 




-0.009 


-0.011 


-0.009 


-0.010 




0.053 


0.059 


0.059 


0.051 




-1.249 


-1.241 


-1.249 


-1.264 




0.256 


0.218 


0.259 


0.266 




-0.036 


-0.033 


-0.034 


-0.037 




0.135 


0.157 


0.145 


0.135 




-1.437 


-1.403 


-1.437 


-1.468 


Cio(a'a') 


0.539 


0.459 


0.546 


0.565 



Table 3: Wilson coefficients of the electroweak penguin operators at typical hadronic scales 
fif, = 4.4 GeV and = 1-3 GeV in units of a in the NDR scheme and HV schemes for 
= Myy, = frtt (see text). 



the NDR scheme. Again, the reduction of the dependence on the renormahzation scale 
of the top mass is remarkable. In contrast to C-j and Cg the NNLO corrections to Cg 
and Cio are substantial in a large range of fit and a "naive" choice = in the NLO 
expressions would, in particular in the case of Cg, totally misrepresent the true value of 
these coefficients. This peculiar behaviour of Cg and Cio can be traced back to the fact 
that Cs{M^) = Cio{M^) = at NLO. 

In Table ^ we show the results for the Wilson coefficients as in Table ^ but this time 
choosing fit = frit. We observe a significant reduction of the NNLO corrections in the 
case of C7 and Cg relative to Table The corrections to Cg and Cio in the NDR scheme 
increase and decrease, respectively. In the case of HV they are smaller than in Table § 
but this time Cio is slightly enhanced. In any case the strong scheme dependence of Cg 
and Cio observed in Table ^ is also evident here. 



7.4 Scheme Dependence of Cg and e' /e 

The strong scheme dependence of Cg at the NNLO level is welcome. In the case of the 
CP- violating ratio e'/e, the operator Qg is by far the most important electroweak penguin 
operator due to its large AI = 3/2 matrix element (Q8(/^x))2 = -Bg {fJ'K){Q8{fJ'K))T'^ , 
where "vac" stands for the vacuum insertion approximation. The scheme dependence of 
{Q8{i^k))2 resides fully in Bf^'^\jjLK)- As the contribution of Qg is the dominant 0{a) 
contribution to e'/e one expects that the product B^^'^\ij,k)Cs{ij.k) is approximately /j^k 
and renormahzation scheme independent with small hk and scheme dependences to be 
canceled by contributions of other operators which mix with Qg under renormahzation. 



This is supported by renormahzation group studies [|14[ which also show that at the NLO 



level B^^'^\ij,k)Cs{ij,k) is only weakly dependent on /x^ for 1 GeV < /i^ < 2 GeV. 
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The situation with the scheme dependence of i^g {fJ^K)C8{f^K) is different. Only by in- 
cluding the NNLO corrections to Csinx) calculated in the present paper B^^'^\fiK)Cs{fJ'K) 
turns out to be almost scheme independent, whereas a substantial scheme dependence is 
observed at NLO. Indeed using the results of Table 3 and -Bg^nv — 1-2 B^i^BR HH find 
for fix = 1-3 GeV 

Bi'/^iMCi^'M f 1.29 NLO, ^^^^^ 



Bii^iMCr^'M I 1-03 NNLO. 



This result can be understood by recalling that at NLO Cg has the formal expansion 
0{a/as) + 0{a). Now the NLO term 0{a) is substantially larger than the leading term 
0{a/as) mainly due to the Z°-penguin diagrams which contribute first at the NLO level. 
In evaluating numerically the product B^^'^^Cg one effectively includes a term 0{aas) 
which originates in the product of the large 0{a) NLO term in Cg and the scheme de- 
pendent 0{as) correction in Bf^'^\ As the 0{aas) term in question is really a part of 
the NNLO contribution and moreover it is substantial, the resulting scheme dependence 
of B\ Cg at NLO is large. Including the 0{aas) corrections to Cg removes this scheme 
dependence to a large extent as seen in (|7.14|) . 

We would like to remark that the corresponding product Sg^^^^Ce related to the domi- 
nant QCD-penguin operator in e'/e, exhibits a much smaller scheme dependence at NLO 
than 5f /^^Cg. In this case the ratio corresponding to ( |7.14D is found to be LOS at NLO. 



This is related dominantly to the fact that the NLO contribution to Cg is relatively small 
compared to the leading term in contrast to the case of Cg as discussed above. 

What is the impact of our results for Cg on e' /e ? Clearly the main theoretical uncer- 
tainties in e' / e reside in the values of the hadronic matrix elements which are substantially 
larger than the renormalization scheme uncertainties just discussed. Yet our calculation 
of NNLO corrections allows us to reduce considerably the /i^ and in particular the renor- 
malization scheme dependence in the electroweak penguin sector. However, in order to 
give a shift in e' /e due to NNLO corrections one would have to include similar corrections 
to QCD-penguin contributions and subdominant 0{aas) terms. 

On the other hand, the inspection of Table 3 and ( |7.14]) shows that the role of the 
electroweak penguins for fixed hadronic matrix elements is increased by roughly 16% in 
the NDR scheme and decreased by roughly 7% in the HV scheme compared to the NLO 
results. As electroweak penguins contribute negatively to e' /e, which is dominated by 
a positive contribution from the QCD penguin operator Qe, the NNLO corrections to 
Cg calculated here suppress (eV^)NDR and enhance (£:V^)hv over their NLO values. As 



an example taking central values of the parameters used in ||38| and including NNLO 
corrections to Ct-io^hk) we find 

1 5.9.10- NDR, 
[ 6.3-10"^ HV. 

to be compared with 7.1 ■ 10"^ (NDR) and 6.1 ■ 10"^ (HV) at NLO. Here in contrast 
to we have used -Bg^nv — l-^-Bg^NOR and -Bg^HV — l-^-Bg^NDR which results in higher 
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HV values than obtained there. We have checked that the remaining scheme dependence 
resides dominantly in the QCD-penguin contributions for which NNLO corrections are 
unknown. For larger (smaller) B^^'^^ at fixed Bq^^'^^ the impact on e'/e coming from 
NNLO corrections to electroweak penguin contributions is larger (smaller). 



7.5 Universality of the Functions X and Y 

As discussed in the Introduction, any decay amplitude can be written as a linear combi- 
nation of mt-dependent functions present in the initial conditions Cii^My^,). In the absence 
of QCD corrections the gauge independent set relevant for non-leptonic and semi-leptonic 
rare K and B decays is given by p[ 



^0 — Cq — ABq , Yq — Co — Bq , Zq 



Cn 



4 



(7.16) 



and Eq, with Cq,Bq, Dq and Eq entering Cj(M^y) in ( p.9|) -( pTT5| ). Here we will only discuss 
Xq and Yq. In the case of semi-leptonic FCNC processes the inclusion of 0{as) corrections 
to Z'^-penguin and box diagrams generalizes Xq and Yq to 



Xe{xt) = Ce{xt) - ABe{xt, +1/2), Y,{xt) = C,{xt) - B,{xt, -1/2), 



(7.17) 



where Ce{xt) = C{xt) is given in ( |5.2| ) and 



a.. 



B,{xt, ±1/2) = Boixt) + -^Biixt, ±1/2) 

Arc 



(7.18) 



with _Bi(xt,±l/2) given in ||, 0. Concentrating first on the operators and Qg and 
0{aw) terms in (|2.9|) and ( p.l5|) , respectively, our calculation of gluonic corrections to 
box and Z°-penguin diagrams provides the generalization of Xq and Yq relevant for non- 
leptonic decays as follows 



ACsiM^) = -^[2Yq-Xq] 



-^[2Yq + 2Xq] 

DTT 



-'^12Y, + 2X,], 

DTT 



where 



X,{xt) = XQ{xt) + ^(Ci{xt)-^G{xt,x,) + B'i{xt)^ 
Y,{xt) = YQ{xt) + ^(c^{xt)-^G{xt,x,)-Bfixt) 



Analogously we can write in the case of the operators and Qiq 



AC4(Mv, 



2Yq — Xq 



ACio(M,; 



2Yq + 2Xq 



where 



a. 



An \2 



1 



^^i^t) = 7^ U ^(^f + Bii^t) , Yq{xt) = 7^ - G{xt, X,) - Bfixt) . (7.24) 



47r V2 



(7.19) 
(7.20) 

(7.21) 
(7.22) 

(7.23) 
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fit = Mw 


fj't = m 




NDR 


HV 


NDR 


HV 




0.912 


0.894 


0.980 


0.962 


Y 
% 


0.911 


0.908 


1.006 


1.003 


Vgi 


0.985 


0.968 


0.986 


0.966 


Y 


0.993 


0.991 


0.994 


0.990 



Table 4: rj factors for the functions X and Y in different schemes and for different /it. 



Evidently, at NNLO in non-leptonic decays more mt-dependent functions appear than in 
the case of semi-leptonic FCNC processes. Moreover additional functions are necessary 
to describe the mj-dependence of the coefficients C7 and Cs as seen in ( [7.3| ) and ( [7.4|) . 
Furthermore, gluon corrections to photon penguins and electroweak corrections to gluon 
penguins will introduce new mj-dependent functions not present in semi-leptonic FCNC 
decays. 

We conclude therefore that at the NNLO level in non-leptonic decays the structure of 
mt-dependence is much more involved than in semi-leptonic FCNC decays. On the other 
hand, if we restrict our discussion to the dominant mt-dependence residing in Cgi^My^,) we 
can say something concrete about the violation of the universality of the mt-dependent 
functions addressed briefly in the Introduction. To this end we write 



Xg{xt)=Xe{xt) + -^ 

4:71 
CV 

Y,{xt)=Y,{xt) + -^ 
47r 



AB^{xt-l/2)-\G{xt,x-,) + B^{xt) 
6 

fii (xt,-l/2) - \G{xt, x-,)-Bi{xt) 
6 



= m^X,(xt)=r7fXo(xt), (7.25) 



= r;5r,(^t) = <l?.(xt). (7.26) 



Clearly, the size of the various 77 factors depends on the choice of /Xf In Table ^ we report 
their values for yUt = and fit = in the NDR and HV schemes. As usual, we fix 
livv = Myy. For fit = M^,- the universality of X and Y is broken at 0{as) by terms 
which are relatively small with respect to the NNLO correction, although not negligible 
in the HV scheme. This follows also from our previous remark that for this choice of 
scale the largest contribution to Cg comes from Ci{xt), which is the same for hadronic 
and semi-leptonic decays. In the case of fit = nit, however, the C(as) corrections to X, Y 
never exceed 4% and Ci{xt) plays no longer a dominant role. Although the universality of 
X and Y is broken by effects which are of the same order of the NNLO correction, these 
corrections are anyway much smaller in this case for X and Y than in the case fit = Myy. 
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8 Summary 



In this paper we have calculated the 0{as) corrections to the Z°-penguin and elec- 
troweak box diagrams relevant for non-leptonic AS" = 1 and AB = 1 decays. This calcu- 
lation provides the complete 0{awO:s) and 0{awC(sS'iia'^ 9wTnf) corrections to the Wilson 
coefficients of the electroweak penguin four quark operators relevant for non-leptonic K- 
and 5-decays. We have given arguments supported by numerical estimates that the cor- 
rections calculated by us constitute by far the dominant part of the next-next-to-leading 
(NNLO) contributions to these coefficients in the renormalization group improved pertur- 
bation theory. 

The main results for 0{as) corrections to the Z^-penguin diagrams can be found in 
{ f).8\) and ( |5.21|) . Those for the box diagrams in ( |6.19| ) and ( |6.22|) . The main results for 



the Wilson coefficients of the electroweak penguin operators are collected in ( [7.3|) -( [n3D . 
The numerical values of these coefficients are collected in Tables 1-3 and in Figs. 5 and 6. 

Our main findings are as follows: 

i) The inclusion of NNLO corrections allows to reduce considerably the uncertainty 
due to the choice of the scale fit in the running top quark mass mt{fit) present in 
NLO calculations. This is illustrated in Figs. 5 and 6. 

ii) While NNLO corrections to Ci and Cg are generally moderate and very small for the 
choice /it = rrit, they are sizable in the case of Cg and Ciq. This is illustrated in Tables 
5 and 6. In particular we observe substantial renormalization scheme dependence in 
Cs and Cio, whereas the scheme dependence in C7 and Cg is significantly smaller. 

iii) The strong scheme dependence of Cg allows to cancel to a large extent the scheme 
dependence of the matrix element {Qs)2 relevant for e'/e so that the contribution 
of this dominant electroweak operator to e'/e is nearly scheme independent. This 
should be contrasted with the existing NLO calculations of e'/e which exhibit sizeable 
scheme dependence in the electroweak penguin sector. 

iv) In the case of AB = 1 decays the most important among the electroweak penguin 
operators is the operator Qg. As the NNLO corrections for Ht = fnt are in the ball 
park of a few percent, our results have smaller impact on non-leptonic AB = 1 
decays except for the reduction of the /^rdependence. 

v) We have also investigated the breakdown of the universality in the mj-dependent 
functions X and Y. As these functions are dominated by the contribution of the 
color singlet Z°-penguin diagram which is universal, the breakdown of universality 
through color non-singlet Z^-contributions and box diagrams is small as illustrated 
in Table 4. 

Although we have seen that there are arguments suggesting that our subset of NNLO 
corrections is dominant, several other contributions have to be calculated in order to 
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complete the NNLO analysis for non-leptonic decays. We have discussed this formally in 
Section 3. A step in this direction has been made recently in |[T^ where 0{a'^) corrections 
to the initial values Ci_6(^m/) have been calculated. Yet the complete 0{as) corrections 
to the photon penguin diagrams relevant for non-leptonic decays and in particular the 
three loop anomalous dimensions 0{al) and O(aa^) of the set Qi-io are unknown. The 
present work and the complementary calculation in |]13| constitute the first steps torwards 
a complete NNLO calculation of non-leptonic decays and we have demonstrated here that 
the NNLO corrections to the Wilson coefficients of electroweak penguin operators are of 
phenomenological relevance. 
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Appendix 



In this Appendix we report the analytic expressions for the functions ZQ^i{xt,Xz) intro- 
duced in ( ^.12| ). They read 



Xt (20 - 20x2 - A57x^ + 19xtx^ + 8x1 



32{xt - l)x^ 

Xt {lOx'f - llx^x-, - xt(30 - IQx-,) + 4(5 - 17x^ + xD) 

iQ{xt - lyxz 

Xt (10 — lOx^ — ITxz — XtXz — 4:xl) 



In Xt 



16{xt - l)xz 



XzilOx^ - xt{4:- Xz) + 8xz) , 2 
mXz TTTT TTT^ In Xt 



32(xi - 1)2 



-'-2 12 

-x^m X, 
4 ' 



+ I2xt + xf 5{xt-iy{2 + xt) 



4x, 



16x2 



(Al) 



12 - 3x? - 3x^2(4 - Xz) + 4xt(3 - x^) + Ax^ - x 



In Xj In x, 

K^t - 1)2 

12xt + x2 5(xt- 1)2(2 -fxt) 3 (4 + 8xt + 2x2 - X?) 



V 2x, 

{xz — 1)2 (5 — 6xz — 5x2) 



8x2 



Li2(l - Xz 



4(Xi-l)2 

5 - 16x^ + 12x2 + 2xt 



4x2 4x2 
xt(4 - Xz) (88 - 30x^ - 25x2 - 2xt(44 - 5x^ - 6x2)) 



j Li2(l - Xt) 
C(2) 



32{xt - l)2x, V 4 

IQxf - xt{20 - Xz)x\ + 8x^ - 8x^^(14 + hxz) + 8xl{\2 - Ixz + xl) / Xz' 

32(xi - l)2x, '^fe 
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22 + 33xt - xl b{xt -l){2 + xt) ^2 + bx] + IQx^ + Xt{l^ + x^) 



16(xt — l)xz 



16x2 



16(xt - 1)^ 



^{Xt,Xz) 



and 



xt(20 - 20x2 - 265a;^ + %lxtx^ + 8x 



48 (xt - l)x, 

xt (lOxf - 15x^23,^ + 4(^5 _ 7^,^ + 2x2) - xj(30 + 20x2 + 4x2)) 

24(xt - l)2x, 
Xt(10 -lOxf -33X;j +15xtX2 -4a;^) 



Inxt 



24(xt - l)x. 



Inx, + 



x^{8 -16xt +2x2 +10x2 + 7xtx^) 



48(xt - 1)2 



X2(4 + 5x2) 2 
H ^^r-; m X2 + 



24 



20 + 6xt + xf 5(xt - 1)2(2 + Xt) 



12x^ 



24x? 



In^ Xt 
(A2) 



3xf + 2x|(12 - x^) - xt(18 - 16x2 + xl) - 2(9 + 4x2 - 

12(xt - 1)2 



In Xt In X2 



/20 + 6xt + x2 5(xt- 1)2(2 + xt) 6 + 6xt-8x2-x3\ 
^ I Q^. 12^1 2(x, - 1)2 j ^''^^ - ""'^ 

(X2 - If (5 - 10x2 - 7x2) 10 - 40x2 + 36x2 + 4x^ + 5x^ 

Li2U - X2j H — -T^ 



6x2 



12x2 



Xt(x2 - 4) (24 - 26x2 - 13x2 - 6xt(4 - X2 - x2)) /x 



+ 



16(xt - 1)2x2 

24x? + 12x2(14 + X2) - 2x2(4 + 5x2) - xt(80 - 36x2 + 7x2 

48(xj - 1)2 

lO-Xt-if 5(xt - l)(2 + xt) 6 + 3x2 + 14x2 + 5xt (7 + X2) 



2x2(2 + xO 
3(xt - l)x2 

X2 
Axt 



l{xt - 1)X2 



24x2 



+ 



24(xt - 1)^ 



$(Xt, X2) . 



We recall that ^(2) = 7r2/6. The function (j)(yz) appearing in the above expressions is given 

by 



^J-^^ Cl2(2arcsinv^) , 



< ^ < 1, 



-4Li2(i^) + 2W(i^)-ln'(4^) + 2C(2) , z>l, 



where Cl2(x) = ImLi2(e'^) is the Clausen function and (5 — ^Jl — l/z. Defining 



A(x, y) = Ml-x-yf - 4xy, 



(A3) 



the function $(x, y) admits two different representations, according to the sign of A2(x, y). 
For > we have 

$(x,y) = A |21n f ^ + ^-^-^ '1 In ( ^^^±1^] _i,^i,y 



-2 Lis 



1 + X — y — A^ 



2Lji^£±l^ +2C(2) , (A4) 
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while for < 



2 arccos 



-1 + X + 



2 arccos 



1 + X — y 



+ 012 



2 arccos 



1 — X + y 



Additional details on this function can be found in 39 



(A5) 
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